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FESHBACH RESONANCES IN POSITRON-HYDROGEN SCATTERING* 

Richard J. Drachman 
Theoretical Studies Group 

Goddard Space Flight Center, Greenbelt, Maryland 20771 

Using a wave function including the ls-2s-2p eigenstates of hydrogen and 
the Is ground state of positronium, we have calculated the position and width of 
the lowest Feshbach resonance in the e + -H system, which lies between the 
positronium rearrangement threshold and the n = 2 hydrogen threshold. The 
method differs from the usual close- coupling method, resembling in its formu- 
lation the projection operator technique, although orthogonality of closed- channel 
and open-channel functions is not demanded. The resonance lies 9.8281 eV above 
the ground state of hydrogen and has a width of 0.029 eV in this approximation. 
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I. Introduction 

In the past few years the theory of non-res onant positron-hydrogen 
scattering has been throughly investigated . 1 For S-waves, the situation below 

2 3 4 5 

the positronium formation threshold is very clear, while recent work ' 
seems to be leading rapidly to a solution of the two-channel problem above 
threshold. 

A number of investigations into possible structure in the S-wave problem 
have also been carried out with negative results. Attempts to find a three-body 
bound state by variational methods 6, 7 were not successful and it was finally 
proven 6 that no such state exists. An indication of a scattering resonance 9 
below the positronium formation threshold is probably an artifact . 7 ’ 10 

The usual method of searching for such compound-atom resonances is 
with the aid of Feshbach projection operators . 11 The idea is to remove the 
open-channel parts of the scattering problem, leaving a residual Hamiltonian 
possessing a discrete spectrum. Resonances of the scattering problem can be 
shown 12 to lie in the vicinity of the eigenvalues of the residual Hamiltonian, 
and in practice such resonances are really extremely close. For this reason, 
it has usually been assumed that resonances exist whenever eigenvalues of the 
residual Hamiltonian appear; the shifts A due to coupling between closed and 
open channels have never yet been large enough to alter the qualitative situation. 
Several attempts to find eigenvalues below the positronium formation threshold 
gave negative results 2,7,14 and it is now believed that no resonances exist in 


this region. 
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It is natural to consider next the possible existence of resonances just below 
the n = 2 threshold in hydrogen. One might be optimistic about finding such 
resonances in the positron -hydrogen case since it is well known that the 2s-2p 
degeneracy produces an infinite sequence of resonances in the electron-hydrogen 
case. 15 Not only do the eigenvalues of the residual Hamiltonian form a clear- 
cut infinite sequence lying wholly below the threshold, 15 but scattering calcu- 
lations 16 imply very small values of A, at least for the low-lying resonances. 

Mittleman 17 has examined the analogous problem for positrons. Here the 
situation differs because a second (positronium) channel is open and it is more 
difficult to handle the non-orthogonal functions representing the two open 
channels. Nevertheless, he showed that the potential is the same asymptotically 
as in the electron case, falling off as the inverse square of the positron co- 
ordinate. Hence, the same conclusion follows: there must be an infinite set of 
eigenvalues here also. No information on the actual energies was obtained, 
however, since that depends on the short-range behavior of the potential as well. 

Also important, however, is the fact that Mittle man’s paper 17 gives no 

information about A, the energy shift, which in this case involves coupling with 

both open channels and might be large enough to move the resonances far away 

from the positions of the eigenvalues. Although this phenomenon has never been 

observed before, it must hot be overlooked as a possibility. This is especially 

true since a fairly close search 5 has failed to find any resonances in a two- 

1 8 

channel lower bound calculation, while a close- coupling calculation omitting 
the positronium channel gave the first three S-wave resonances explicitly. 
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It is the purpose of this Article to clarify the situation with regard to the 
two-channel resonances. In Section n we formulate an approximate two-channel 
scattering equation, from which the lowest resonance position and width can be 
obtained by a modification of the Feshbach operator formalism. In section III 
we discuss an explicit form for the correlation function describing the resonant 
state. Section IV gives our conclusions and a discussion, while some mathematical 
details are in an Appendix. 

II. Formulation of the Approximation 

Let us assume for the moment that some normalized function $ exists and 
describes well the long-lived three-particle state responsible for the resonance. 
(Later on a particular form for $ will be selected). A two-channel scattering 
function can be constructed in the following form: 

= * (*) F(R) 4 (?) + c. § . ( 1 ) 

Here the ground state of hydrogen is designated Is and that of positronium is Is, 
x and r are the positron and electron coordinates respectively, while the two 
coordinates 

g J t~- S-& (U) 

describe the center of mass position and relative coordinate of the positronium. 
The two functions v and F and the constant C are to be determined; the last is 
large when the energy is near resonance. 
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Coupled equations for these three unknowns are derived by standard vari- 
ational methods. The resulting set is as follows: 


( H 


Here 


5 [H-E ]f-=0. 

H = - - ^ + 2 (lif " 7 l ~ t ) 


(2a) 


(2b) 


(2c) 


and 


E = - Eg. + *£•) 


From Eq. (2c) we can obtain an expression for C in terms of ■)( and F as 

c - (£-e)'[$<kk. $ (m-£)<& 4 

_ -i (3) 

+jpf Jr $ (H-E)4> a r] j 

where we have defined 

£ = 5 H $-> <4) 
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and have re-written the volume element in the second integral. The energy £ 
should lie near ihe resonance represented by $ . It is convenient to define two 
new functions 

S (*) = ^4 («-£) <f < S| 

and 

T(R)= jdf $, 

and to rewrite Eq. (3) (after integration by parts) as 

C= (£-£)“' [I^ s + P£ FT], < 

Inserting this expression into Eqs. (2a, 2b) we derive the coupled equations in 
the optical potential form 


■ + <y^>(F-£)'y=a 

where and V are column vectors defined as 



(V • ^ is the quantity in square brackets in Eq. (6), and £ is the 2x2 
matrix whose elements are as follows: 
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H \ (M - i,2) 

where 

2.1 ~ d > Ji - Pj 01 ^ 04S j 02 " 0t$ * 


(The diagonal elements of £ are no more complex than differential operators, 
while the two off-diagonal elements are integral operators understood to operate 
on F and )( respectively. Details are given in the Appendix.) Eq. (7) is thus a 
set of coupled homogeneous integro- differential equations, whose possible 
resonant behavior is due to the second (optical potential) term. The similarity 
to the usual optical potential methods is clear. 

1 9 

Following Chen and Chung we solve Eq. (7) by a non-iterative method. 
Letting 





%*+ <VY> (£-£)"' t 


( 10 ) 


where j = 1, 2 refers to two regular, linearly independent solutions of Eq. (7), 
we must solve 



(Ha) 


£ (f’i + v = o. 


(Hb) 
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Only one particular solution of the inhomogeneous equation (lib) is needed. 
Using the definitions given in Eqs. (8) and (10) one finds 

where the bracket notation is an obvious extension of that defined in Eq. (8). 
Thus from two solutions of the homogeneous equation (11a) one can construct 
the "coupled static" approximation, while through Eqs. (10) and (12) the solution 
of the inhomogeneous equation (11b) leads to two solutions incorporating the ef- 
fect of coupling to the "resonant" state <s>. 

The most convenient way to represent the scattering parameters, for 
coupled and uncoupled cases alike, is in terms of the K-matrix. For any two 

A i 

regular, linearly independent solutions (either *P or H'5 ) one obtains the 
following asymptotic form 



For i, j = 1, 2 the quantities A , and B form 2x2 matrices and the K-matrix 

ij n 

is defined as 



( 14 ) 
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and is related to the S-matrlx by 



(1 +XK)(l-x-K)* 


(15) 


Let us now choose a special set of solutions to the homogeneous equation 
(11a) with the matrix A equal to the identity matrix and hence B = K°, the 
coupled static K-matrix. Noting the fact that the off-diagonal elements of K are 
very small compared with the diagonal elements 3 ’ 4 ' 5 [see Table 1 ], we 
simplify the analysis by setting K° = K° i =0. This approximation is equivalent 
to treating the two physical channels (e + + H and Ps + p) as eigenchannels in the 
coupled static approximation. It does not imply that the coupling between these 
two channels is unimportant, and in our numerical solutions we have not 
neglected the off-diagonal elements of £. (See Sec. IV.) Only in the asymptotic 
region are the two physical channels assumed to be uncoupled. 

Wiih the above simplification it will be convenient to choose the asymptotic 
forms of and f as follows: 



(16a) 


(16b) 


(16c) 
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From the expressions in Eqs. (10 - 14) we can now construct the K- matrix in- 
cluding the coupling to the resonant state as follows: 



To relate the various quantities appearing in Eq. (17) to the resonance 
position and width, we note that the eigenphase sum exhibits typical Breit- 
wigner behaviour. That is, 


V % = H ► C + ^ J (18> 

where 17 and vf are the resonant and non-resonant (static exchange) eigen- 
phases, T is die width and E is the position of the resonance. (The individual 
eigenphases may or may not exhibit such behavior). Using trigonometric 
identities we can show that 

jTq " /i- <x° 2 J [ J 1 - cf * (6jr *0 
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where 


P= T/> K, PV 17, K C ; Qb det Q° = det K° C 

and we recall that the eigenvalues of K are tan 17 2 and tan 77 2 . Using Eq. (17) 
we find that 



( 21 ) 


Inserting these expressions for P and Q into the left side of Eq. (19), we can 
make the following identification 


-p.il + 

4? - 

r 2z - 

r 


(22) 

and hence (from 

Eq. (17)) 




Er- 

ft 


£ + A 


r = 

n 

-L. 

i Ci <i- 

#>. 

(23) 

\ 


l" 1 



These relations 

are the principal results of 

this section of the paper. 

and are 


correct provided (as usual) that 17 ?, E„ and F can be treated as constant across 

J K 
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the resonance, and that A is not too large. We have essentially re-derived the 
results of Chen and Chung 19 for the two- channel special case being examined 

i 

here with special emphasis on the asymptotic forms of JP 0 and ^ In addition, 

we have not required that the resonant part of the trial function be orthogonal to 

20 

the rest of the function as is usually done. 

in. The Close- Coupling Correlation Function 

Since we know that the degeneracy between the 2s and 2p levels of hydrogen 
produces the infinite series of Feshbach resonances 15 ' 17 , it is natural to choose 
the resonant term in our trial function to include those two states: 

$ = N x 1 + nr(-x)(p 2f p) (24) 

where 

(p 2S (A = (32*^ e J (2--^ 

and 


The normalization constant is given by 

z 

(25) 


N - F 4 TT C (uS'v 3 ) 

■^C 


By minimizing £ with respect to u and v in Eq. (4) we obtain the equations 



14 



The potentials are easily obtained: 21 



Assuming a value for the eigenvalue e, the two coupled second-order differential 
equations (26) were solved numerically: two regular, linearly independent 
solutions were obtained by integration outward from the origin and two others 
by integration inward from some large radius x q . Admissable solutions are 
those for which M and dM/dx are continuous at the matching point x; the 

-V ~ 

continuity condition is equivalent to the determinantal condition 



15 


det Q-jl'i o . 


(28) 


The index i varies from 1 to 4 to specify the four linearly independent solutions, 
while j varies from 1 to 4 to denote a = u, v, du/dx and dv/dx, respectively. We 
were able to obtain the first two eigenvalues accurately (with x 0 - 100) using 
Eq. (28). They are 


r Ry 3 

___ Q j. Qf 10 - 


(29) 


In the next section we will discuss these results, comparing them with previous 

« Q 

work. For the moment, let us just note that the ratio of the two eigenvalues 
e j/ € l = 0.0613; this is to be compared with the ratio 0.0583 based on the x -2 
potential alone. 15 

We can now proceed to the calculation of V [Eqs. (5) and (8)] , required for 
obtaining the resonance parameters. Its two components are [from Eqs. (5) 
and (24)] 
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T(R)= iN^d? <^(t) x l [ U(z )[(£- W 

<30b) 

+ 'i r (x)| - _ f- vCj^ / -/i. (p 2S 

Equation (30a) is simple because of orthogonality while Eq. (30b) has been 
simplified somewhat by removing any derivatives that appear through the use of 
Eq. (26). We can re-write Eq. (30a) as 

5 (x) — hi x 1 [V s C x ) a ( x ) + Vf, (x) /vex)"] ^ (S1) 

where 



while Eq. (30b) requires double numerical integration. In Fig. 1 we plot 
S = 2 77 1/2 x S and t = 4 v'V RT to give an idea of their range and form; only the 
first resonance, corresponding to e lf will be discussed from now on. 

IV. Results and Discussion 

Once the coupling term V [Eqs. (8) and (30)] and the operator C [Eq. (9) and 
Appendix] are in hand we need only to solve the two sets of coupled equations 
(11), as discussed above, to obtain the position and width of the lowest e + - H 
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resonance. These parameters constitute the main results of this work and will 
be discussed below. 

In Table I we list the non-resonant scattering parameters, obtained from 
the homogeneous equations (11a), and compare them with the results of Ref. 5. 
These depend on the accuracy of our integration technique and our representation 
of the integral operator £, as discussed in the Appendix. We used a simple, 
less accurate method than that in Ref, 5, truncating the coupled equations at 
x = 20 and R = 12, and it is reassuring that only small errors were thus intro- 
duced. The two off-diagonal elements of K° should be equal; our truncation re- 
sults in significant inequality, although the average of the off-diagonal terms is 
in almost exact agreement with Ref. 5. Note, however, that these off-diagonal 
elements are small, as we have assumed; the two eigentangents are almost in- 
distinguishable from the corresponding diagonal elements of K°. We also show 
results of the uncoupled approximation resulting from the omission of the off- 
diagonal parts of C. It is interesting to note that this omission hardly affects 
the e + - H channel, 4 while the Ps - P channel is modified radically. 

ha Table n we give the properties of the lowest resonance. (Details of the 
numerical calculation are in the Appendix.) The main conclusion is that the 
resonance is quite far below the n = 2 threshold (0,3763 eV) and narrow enough 
to be well defined, with a width only 8% of the distance to threshold. Several 


other points should be noted: 
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1 . When we drop the off-diagonal parts of C the widths and shifts of the two 
channels separate, and one in effect is calculating with the two distinct wave 
functions (Is, 2s, 2p) and (Is, 2s, 2p). Since the decay of the resonant state 

(2 s, 2p) can occur through either open channel the total width and shift are 
sums over the two distinct cases. The first of these (e + - H) should correspond 
to the case treated by Seiler et al. 18 Comparing the first two lines in Table n 
we see that the resonance positions are in quite good agreement, while the 
width is not in agreement at all. We do not understand this discrepancy, but 
have noticed that the ratio of the two results is close to the numerical value of 
the rydberg in eV. It may be that an error was made in Ref. 18 in expressing 
T in eV. Assuming that, we divide their value for T by and obtain a value of 
8.305 x 10 _s eV, differing from our result by only 1.3%. 

2. The second (Ps - P) channel has a much larger effect on both width and 
shift than the first. In the uncoupled approximation this amounts to a factor of 
about 7 in A and about 35 in r. This is expected for A, since only the sum 

£.+ A has physical meaning and £ itself depends on the separation between "open” 
and "closed" parts of the wave function. In channel 1 we have used an open- 
channel function orthogonal to the (2s, 2p) closed- channel function, but not in 
channel 2. We expect E R = £+ A to be fairly reliable and probably an upper 
bound to the true resonance energy. On the other hand, T is an observable it- 
self, and its sensitivity to the several approximations shown perhaps implies 
that even our best result may be unreliable: our best r is about 10 times 



larger than our uncoupled width, while our best A i s only 1 .4 times smaller 
than our uncoupled shift. 

3. In order to compare our non-resonant results with previous work, 5 we 
solved the scattering equations at k = 0.85 instead ofk = 0.847781 which cor re 
sponds to the exact value of £. This is certainly an insignificant error. An 
improvement would be to calculate A again at the value of k corresponding to 
E r , continuing to recompute A until convergence is attained. This correction 
is also expected to be small. In fact, the value at our E„ is k = 0.8499; by 
accident we have evaluated A and r at very nearly the best possible k. 

4. The resonance described here may be compared with the lowest 
resonance in e - H scattering, as treated in the ls-2s-2p approximation 18 as 
follows: 

— J (A Re _ k' 2 %2p S ) = O- icO 



o. 5 H . 


Thus, the two two types of resonances are not qualitatively different. The long- 
range potentials are the same and the positronium channel is similar to the 
exchange channel, while the change in sign of the short-range static potential 
raises the e + - H energy and lowers the e - H energy. 



APPENDIX 


A. Evaluation of the operator C 


In order to solve Eq. (7) [or its equivalent, Eq. (11)1 we must evaluate the 
integro-differential operator P, [Eq. (9)]. Its diagonal elements are very easy 
to obtain: 




J. 




-A + ^is ( . x )j 


(Al) 


^21 - Pis 




£ 


-%-H] <P t 


<■?) 


(A 2) 




Here we have used either (x, r) or (R, p) variables as convenient and have made 
the definitions 


li = -'t- r ; H, (?) r -2 7/ - i (A3) 

also note that <p ls and are eigenfunctions of these Hamiltonians with eigen- 
values E 1S = -1 and E js-=“ 1/2 respectively. The potential in Eq. (Al) is 
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the static or Hartree potential felt by the positron in the field of an undistorted 
hydrogen atom. (The corresponding potential in the positronium channel vanishes 
due to antisymmetry and is absent from Eq. (A2). That is why K° 2 = 0 in the 
uncoupled approximation in Table II). 

The off-diagonal elements are discussed in Chan's thesis [Ref. 5] and will 
be briefly reviewed here. The (21) element can be written as 



where we have used the fact that is an eigenfunction of H 0 (r) and have com- 
muted it to the left side. By changing variables from R,p to R, x and integrating 

'Nj 'V 

by parts we obtain the form 




+ x )S'2r\ 


-I s ] 


(A6) 


= Jci% G(S, *)*(*) 

where the factor of 8 is the Jacobian of the transformation. In analogous fashion 
we can write 

F - + 5 T ~ \f^]\ 


- ^8 & (^ 8 ) 


(A7) 
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We can show that G (x, R) = G (R»x) by changing variables from x, R, to r, P 
and writing 

&(&£)- & ( 8>Z) ~ fo(*)-£is - H fl (?) +• E& ] <f> ±& ($ cj>- s (§) r= O, (A8) 

Then we can evaluate the kernel itself most conveniently as 

& -r * & )= +[' 3 ^- s -A- r 

*x “a" 5 1 (A9> 

* +£ *±+$- 

The cross-term is done explicitly and we use the relations between the 
energies tEq. (2) ' . Finally, we obtain the result 

f~ -A - § - 1 "M*] &is ("0 (f). < A10 > 

where r = 2R - x and p = 2 (R - x) . 

Since we are interested only in total S-wave scattering, we allow X and F 
to depend only on the magnitudes x and R. This allows us to carry out the 
angular integrations in Eqs. (A6) and (A7) over jj. = X • R by numerical Gaussian 
quadrature. The final expression of Eqs. (11) involving the operator £ is then 
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=M s -f 4- JdRK(R )Xn(R) + [o,5 (flt) ] 

§ + %\ = 2 K ( R A)f (*) + [«. 

where 

$(x)= , CJ(R)= F(r)/r^t j 


and S, T are defined following Eq. (32) and shown in Fig. 1. The S-wave kernel 
is 


K ( R,x) = S? JT R x ^ [ A’ f f 8 J ‘ + x 1 1 - C A" ‘ 6 ' - 1 - jit ] £ 


-(*♦&) 


(A12) 


where 

o/a 

A = (4R*+x a -4R^), 

0 = + 

C ■= 2R 1 +x 2 -3Rayt'. 

It is easily shown that for small x, K(R,x)cce _3E while for small R, K(R, x)cce" 2x: . 
From these results it appears that K extends further in x than in R. A trouble- 
some ridge occurs for x « R, along which K decreases very slowly. Following 
Chan 5 we have evaluated K(R, x) separately in the four regions bounded by the 
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straight lines x = 0, x = R, x = 2 1/2 R, x = 2R and R = 0, using simple trans- 
formations to avoid singularities and employing numerical Gaussian quadrature. 

We chose to solve the integro-differential equations by converting them to 
a set of linear algebraic equations which is numerically inverted. Specifically, 
we used a three-point central-difference formula to approximate the second 
derivative and the trapezoid rule for integration. A considerable improvement 
in accuracy is obtained by a simple modification; for example, the first homo- 
geneous equation of (All) becomes 



where 


^ 5 (Q, V* „ 



and h is the appropriate lattice interval Ax or A R. We obtained sufficient ac- 
curacy for the homogeneous equations using the following set of intervals: 
AXcQi ft) AX-O.2 ^ AX = 0.4 

\rfi o; ARr =0.06 ^ AR -C.iZ 

h AR = 0.24 fa 
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This gives a total of 100 values each for f n and and a set of 200 linear equa- 
tions to invert. Beyond this region, the equations (All) effectively uncouple 
since K ® 0, but the inhomogeneous terms are not negligible. The inhomogeneous 
equations were continued out to x, R = 50 as ordinary differential equations, and 
all the integrals needed to obtain A and r were also carried out that far. 

On the other hand, Chan and Fraser 5 chose to solve the homogeneous Eqs. 
(All) by first converting them to integral equations. This involves some addi- 
tion work, but enables one to use Gauss-Laguerre integration which can repre- 
sent the unknown functions on an optimal lattice of points. Judging from their 
results and ours, it would seem that the integral equation method has some real 
advantages in accuracy and economy. 
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Table II. Properties of the lowest resonance. Our results are in eV, 
using R ro = 13.60583 eV [Ref. 22]. (The results in Ref. 18 were re- 
ported in eV; but the value of the rydberg used was not stated). The 
unshifted energy £ = 9.77896 eV. Our best result is given on the last 
line. 



Width r 

Shift A 

Resonance 
position Ejj = £ +A 

e r - H channel only 

8.20 

(-5) 

8.94 (-3) 

9.7879 

Ref. 18. 

1.130 

(-3) 

- 

9.7845 

Both channels (uncoupled) 

3,03 

(-3) 

6.88 (-2) 

9.8478 

Both channels (coupled) 

2.94 

(-2) 

4.91 (-2) 

9.8281 

1 
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— " ~~ 1/2 

Figure 1. The resonant coupling functions S = 2tt x S(x) and T = 4-n RT(R), 
plotted against x and R respectively. 






